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Abstract: When the D(p + 4)-brane (p = 0,1,2) with delocalized Dp charges is put 
into equilibrium with a spherical thermal cavity, the two kinds of charges can be put into 
canonical or grand canonical ensemble independently by setting different conditions at the 
boundary. Using the thermal stability condition, we discuss the phase structures of various 
ensembles of this system formed in this way and find out the situations that the black 
brane could be the final stable phase in these ensembles. In particular, van der Waals-like 
phase transitions can happen when DO and D4 charges are in different kinds of ensembles. 
Furthermore, our results indicate that the D(p + 4)-branes and the delocalized Dp-branes 
are equipotent. 
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1 Introduction 

Branes are non-perturbative objects in String theory and study on them is helpful in under¬ 

standing the non-perturbative properties of String theory. Various brane configurations are 

also useful in the applications of gauge-string duality. Under the assumption of AdS/CFT, 

study on the near horizon geometry of D3 brane, i.e. AdSs x S' 1 * * * 5 , reveals the properties of 
corresponding strong coupling conformal field theory in the large N c limit at the boundary 
of AdS, 5 , which can also be extended to the asymptotically AdSs spaces. A typical example 
is that the well-known Hawking-Page phase transition of the AdS black hole corresponds 
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to the confinement-deconfinement phase transition in SU(N C ) gauge theory at large N c [1], 
which aroused the recent interests on the black hole thermodynamics in the asymptotically 
AdS space. 

As the solutions of the supergravity, black branes, like black holes, can have their own 
thermodynamics. Study on the thermodynamical phase structure of the black branes is also 
valuable in understanding the non-perturbative nature of String theory. Since black branes 
are asymptotically flat, they are unstable by themselves due to their negative specific heat 
capacity and the Hawking radiation. To study their phase structures, one can put them into 
a spherical cavity which is considered as a reservoir to form a thermal equilibrium, following 
the approach of York in discussing the thermodynamics of black holes [2-4], For different 
boundary conditions at the boundary of the cavity, there can be different ensembles for 
charged black branes. Along these lines, in [5], the phase structure of black p-branes in the 
canonical ensemble was studied in ZD-dimensional space-time, where the temperature, the 
volume of the brane as well as the cavity, and the charges of the brane are fixed. There can 
be van der Waals-like phase transition in this system for d > 2 (p < 5) where D = d + d + 2 
and p = d— 1: a first order phase transition between a large black brane and a small one can 
be found for charge q < q c , while for q > q c there is only one black brane phase, with q c being 
the critical charge at which there is a second order phase transition. This is similar to the 
van der Waals-like phase transition found in charged AdS black holes [6, 7], asymptotically 
flat black holes as well as dS black holes [8, 9]. In contrast, for d < 2 (p = 5,6), there 
is no van der Waals-like phase transition. For uncharged black branes, Hawking-Page like 
transitions between black branes and the “hot flat space” can also happen in the canonical 
ensemble similar to the uncharged black hole case [2, 3]. The phase structure for black 
branes in grand canonical ensemble is different from that in canonical ensemble for d > 2, 
where the potential is fixed while the charges are not. There is no van der Waals-like phase 
transition, but the Hawking-Page-like transitions between the black branes and the “hot 
flat space” can happen. Bubble solutions[10, 11] which can be obtained from black branes 
by double wick rotations also play a role in the phase structure of black branes, since they 
have the same boundary condition as the black branes. In [12] and [13], bubbles were found 
to be the phases of black branes in canonical ensembles and grand canonical ensembles, and 
can have phase transitions with black branes, and thus enriches the phase structure of black 
branes. In all these discussions, there are no van der Waals-like phase transitions in grand 
canonical ensembles, whereas they do exist in canonical ensembles. The absence of van 
der Waals-like phase transitions in grand canonical ensemble is also true for charged black 
holes either in flat space[4, 8] or in AdS space [6, 7], but they may exist for Gauss-Bonnet 
AdS black holes[14-16] for certain Gauss-Bonnet coupling constants. 

It is well-known that branes can be combined to form composite states or intersect¬ 
ing branes. Among others, Dp-D(p + 4) bound states as solutions of super gravities can 
be constructed by smearing D p branes inside D(p + 4) branes. In extremal cases[17-19] 
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which obey the harmonic function rules, the solutions preserve 1/4 of the supersymme¬ 
tries. There can also be more generic non-supersynnnetric Dp-D(p + 4) solutions in ten¬ 
dimensional supergravities [20-22], These brane solutions are useful in gravity-gauge duality 
applications [23-28] and in constructing lower dimensional black holes [21]. The dynamical 
stability of Dp-D(p + 4) was discussed in[29] . We are interested in the thermodynamical 
phase structure of these solutions in this paper. Since these solutions involve two kinds of 
charges coupled with different RR potentials, the phase structure is expected to be richer. 
One can form different ensembles by fixing either the charge inside the cavity or the po¬ 
tential at the boundary for Dp or/and D(p + 4) branes , and therefore Dp and D(p + 4) 
can be in canonical ensemble or in grand canonical ensemble independently. In [30], the 
phase structure of the D1-D5 system in which both branes are in the canonical ensemble 
is discussed. It has been shown that the smeared Dl-brane alone shares the same phase 
structure as the D5-brane for which the van der Waals-like phase structure can not be 
found. However the D1-D5 combined system displays a van der Waals-like phase structure 
for some region of the charge combinations. And the phase structure is also symmetric 
under the exchange of the D1 charges and D5 charges. One may also be curious about what 
will happen when the two kinds of charges are in different ensembles. In present paper we 
will study the phase structures of the other three combinations of the ensembles: both Dp 
and D(p + 4) in grand canonical ensemble, Dp in canonical and D(p + 4) in grand canonical 
ensemble, and Dp in grand canonical and D(p + 4) in canonical ensemble. We will find that 
the phase structures are also symmetric when one swaps the boundary conditions imposed 
on Dp- and D(p + 4)-branes. For mixed ensembles with one charge in canonical ensemble 
and the other in grand canonical ensemble, one may expect in some charge combinations 
there could be van der Waals-like phase transitions like the canonical ensemble and in some 
other regions there are no van der Waals-like phase transitions as in grand canonical en¬ 
semble. This indeed happens in the D0-D4 system, which has the richest phase structure. 
In D2-D6 and D1-D5 systems, unlike in the D1-D5 canonical ensemble, in the other three 
ensembles there are still no van der Waals-like phase transitions. In our discussion, we only 
consider the stability condition of the brane under a small change of the horizon size, which 
is equivalent to the positive specific heat capacity condition. Nonetheless, we will briefly 
comment on the effect of the more generic stability conditions — the electrical stability 
conditions for both charges — on the phase structure. In the meantime, the technical 
details on this electrical stability conditions is arranged in another paper [31]. 

The paper is organized as follows. In section 2, we review the Dp-D(p + 4) solution and 
evaluate the classical Euclidean actions or thermodynamic potentials for different ensem¬ 
bles. In section 3, we discuss the phase structures of the Dp-D(p + 4) system in different 
ensembles. In section 4, we give a brief discussion on more general thermodynamic stability 
conditions. Section 5 is the conclusion. We also gather some detailed calculations in the 
Appendices. 
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2 The Dp-D(p + 4) brane system 
2.1 The action 

We consider a gravitational system bounded by a big spherical reservoir which can be 
regarded as a spherical boundary at the transverse radius p \,. At the center of the system 
is a pile of parallel Dp-D(p + 4)-branes. The total Euclidean 10-dimensional supergravity 
action in Einstein frame can be expressed as a sum of several contributions [5], 


/ — /EH + I(f> + Ip + /p+4 T /boundary 


( 2 . 1 ) 


where the first term is the usual Einstein-Hilbert action, 

/eh = — J d 10 Xy/gR, ( 2 . 2 ) 

the second term is the contribution from the dilaton field, 

h = ^2 I d 10 x^gd^d^, (2.3) 

the third and fourth terms come from the Dp- and D(p + 4)-brane form field action, 


Ip ~ 4 k 2 
Ip+4 = 4^ 


d w x^/g 

d W Xyfg 


B p2 

(p + 2)! Ip+ 2 ]’ 

(p + 6)! F[p+6] ’ 


(2.4) 


where a n = an d the last term /boundary could admit several boundary integration 

terms, which depends on what ensemble we are interested in. We will come back to this 
term soon. In the above formulae, the constant coefficient n is defined as k = y / 87rG'] ! 0 in 
which Gio is the 10-dimensional Newton’s constant; the vacuum expectation value of the 
dilaton field 4>(r —>• oo) has been shifted to zero. The integration is performed within the 
boundary and outside the horizon if there exists one. 

Now we deal with the boundary term in equation (2.1). There are three terms that 
may contribute to the boundary action /boundary, 


/gh = 

/fo,p — 


/fe,p+4 


^2 j> d 9 x^ (K - K 0 ), 

1 f n g ffi p+40 


(2.5) 


The first action above is the Gibbons-Hawking surface term [32], in which K is the trace 
of the extrinsic curvature defined as 


Rfll/ - 2 T V[y7T^) 


( 2 . 6 ) 
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where n p is the normalized space-like vector normal to the boundary. Kq i n that term 
is defined in the same manner as K but with the metric replaced with flat metric. This 
subtraction term is included to make /qh vanish for flat metric [32] . The second and third 
actions in (2.5) are boundary contributions from the Dp- and D(p + 4)-branes respectively. 
7 in these equations is the determinant of the induced metric on the (4 — p)-dimensional 
boundary, and the (p + l)-form and (p + 5)-form fields A\^ + i^ and A \, p _|_ 5 j are the Rarnond- 
Ramond potentials of Dp- and D (p + 4)-branes respectively. 

The Gibbons-Hawking term should always be included in /boundary regardless of what 
ensemble we are talking about. However, I^ p and h , P +4 are supposed to cancel additional 
boundary terms when doing variations with respect to the gauge field potentials. If we fix 
the gauge field strength F on the boundary and the potential A are not fixed, after partial 
integration the additional boundary terms would emerge, thus we need to include /^ p or 
h,p +4 whose variations would cancel these terms. On the other hand, if we fix the gauge 
potential A on the boundary, there will be no additional variation terms, therefore, 
or Ib,p+A will not be needed. We summarize the above analyses in Table 1. The relation 


Dp 

D (p + 4) 

Boundary action /boundary 

Fixing A [p+1] 

Fixing A [p+5] 

/gh 

Fixing F [p+2 ] 

Fixing A [p+5] 

/gh + Ib,p 

Fixing A [p+1] 

Fixing F| p+6] 

/gh + h,p+ 4 

Fixing F [p+2] 

Fixing F [p+6] 

/gh + h, P + h,p+ 4 


Table 1 . The relationship between boundary conditions and the boundary action 


between boundary conditions and ensembles will be addressed in next subsection. 

2.2 Black brane solution 

The generic non-supersymmetric Dp-D(p + 4) brane solution [20-22, 29] can be obtained 
by directly solving the supergravity equations of motion or by a series of duality and boost 
operations from D(p + 4)-branes. We use following coordinates to describe the solution, 

,Xp , ■■■ , Xp+4, p, (f>i, ■ ■ ■ , (f>4—p) ( 2 - 7 ) 
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where p is the radius in the transverse directions perpendicular to x m , m = 1, • • • ,p + 4. 
In this coordinates the solution reads [30], 


ds 2 


A \p+ 1] 
F [p+ 2] 
A \p+ 5] 

F \p+ 6 ] 

where 


1 —p p—7 


p +4 


A_ 4 A* 8 I A + dt 2 + A_ ^dx 2 + A* ^ dx 


+A 


V — 1 n-4-1 / -7 „2 

4i3^y A s ( d P 


i =1 
„ 2 jo 2 


j=p+l 


A + A_ 


+ p dll 4 _ , 


p—1 3—p 

A„ 2 A* 4 , 


— i- 


A„ - A_ 
A* l A* - A i 


1/2 


dt A dxi A • • • A dx„, 


3-pl / A + \ 1/2 / A_\ 1/2 , , , 

i-— 1--— 1--— dt A dp A dx\ A • • • A dx„, 

p A* V A* / V A* J 

(1 - A \ ^ 2 

— /A + I - -— 1 dt A dxi A • • • A dx p+ 4 , 

\i-A + y 

i - (1 — A + )V 2 (1 — A_) 1//2 dt A dp A dxi A • • • A dx„_i_ 4 , 

P 


( 2 . 8 ) 


A ±(d) 

A*(p) 


1 

1 


3—p 

p± 

«3-p ! 


q 3—p ’ 


P+> P-> 0, 
p i _ p >k > — oo. 


(2.9) 


Here, p can be 0,1,2 for this supergravity solution to describe the Dp-D(p + 4) system. The 
constants p + and p_ are the coordinates for the outer event horizon and an inner event 
horizon respectively, the latter being a curvature singularity. Therefore, the requirement 
p + > p_ is to avoid a naked curvature singularity. The other constant k can be either 
positive or negative, and when it is positive, there is another inner event horizon at p = 
A; 1 /(3—p) which is a curvature singularity as well. When k = p'l p , we have A* = A_, 
and this corresponds to the brane configuration where all Dp-branes are removed from 
the system while D(p + 4)-branes are retained. Later we will see in (2.11) that we require 
p^T P > k in order to make the Ramond-Ramond charge real. 

The Dp - and D(p+4)-brane charge densities can be obtained by integrating the Hodge 
dual of their Ramond-Ramond field strengths, 


Qp 
Qp +4 


l L e “^* F ip+2]0&), 

V 'ZK J Xn 4 -'\ , • ,CCr)4-4 J S' 4 P 


V2kj 

g“p+40 * pjp +6 ](p 6 ), 


— I 




K J S 4 ~P 


( 2 . 10 ) 


where 0 = <t>(pb) and the — % factors come from Euclideanization. The charge densities can 


- 6 - 












be explicitly expressed as follows, 




Qv — 

Sp+4 = ( 3 -P)»4-, eW / 2 ( p + p_)V, 


V 2 i 




( 2 . 11 ) 


where 


and 


^4 — \J9x p+ \x v+ \ (Pb) ' ' ' ffi p + 4 *p +4 (P&) ^4 > ^4 — J dx p+ \ • • • dx p4 - 4 , 


in which 


-3 —p 

A± = A ±(p b ) = 1 - 

Pb 

k 

A* = A *{pb) = 1 — _3_„ j 
Pb 


- v , ~f _ P+ 1 P+1 

P+,-,6 = p +) _, fc A5^ A* 16 = p + ,_, b e^' 


( 2 . 12 ) 


(2.13) 


p+i; 


are the physical radii of horizons and the boundary, and k = ke 4 T 

To make these relations more concise, we define two reduced charge densities, 


Q P = 


y/2.K,Q p e- a P& 2 
(3 - p)V/&i- p pl~ p 


A i 


A* 


1/2 


A_ 


= 1 -^ 1 -^- < 1 , 


A* 


1/2 


Qp+4 — 


'V2 K Q p+4 e- a P+^ 2 \ 3 ~ p . _ 

(3- P )ft,_p J =' / ^ F - 


(2.14) 


With (2.14) we can express p- and A* through p + and these two reduced quantities, 

Ql +4 


p- = 


A* = 


p+ 


A+ + A_ ± a /(A_ - A+) 2 + 4Q2 A+ A_ 

2(1 - Ql) 


(2.15) 


When k = p'l p . the Dp-brane charge vanishes and Q p = 0, so 

2 

Therefore, we should choose in the above relation the “+” sign only, i.e., 


A* = 


A+ + A_ + a /(A_ - A + ) 2 + 4Q2 A+ A_ 

2(1 - Ql) ' 


(2.16) 
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Dp-brane 

D(p + 4)-brane 

Boundary action /boundary 

Canonical 

Canonical 

/gh + h, P + h , P +4 

Grand Canonical 

Canonical 

/gh + h,p+4 

Canonical 

Grand Canonical 

/gh + h,p 

Grand Canonical 

Grand Canonical 

/gh 


Table 2. The relationship between ensembles and the boundary action 


From (2.11), we can see that fixing F at the boundary is to fix the charge density Q 
within the boundary. It is well-known that fixing the charge of a system while allowing 
radiation means we are considering the system in canonical ensemble. If we fix the gauge 
potential A instead, the field strength at the boundary will be changeable, and we will be 
talking about grand canonical ensemble. Now that we have two kinds of gauge potentials 
and field strengths, it is reasonable to assume there are two kinds of ensembles for each 
field, and thus we have the following table (Table 2). 


2.3 Temperature and conjugate potentials 

Like black holes, due to Hawking radiation [33], black brane has a temperature, Tjj, which 
can be easily calculated by the requirement that the Euclideanized metric (2.8) has no 
conical singularity. Then the Euclideanized time t is cyclic with a particular period 

iTTpb (l-A + ) 1 /2(A*-A + ) 1 /2 


(3* = 


2 -y 

(A_ - A + ) 3 -p 


(2.17) 


3 — p 

This period is just the inverse of the Hawking temperature observed by an observer at 


P = oo, 

T h = 1/P*- 

For a local observer at p = p&, the local temperature would be 

f = /3*a|aT A? = |^(1 - A + b (§±) “ (l - 


P-2 

3 -P 


1-A 

A* 


(2.18) 


^ ■ ( 2 - 19 ) 


In order to study the grand canonical ensembles, we need to define two potentials 
which are conjugate to the brane charge densities. So we define to be the potential 
conjugate to the Dp charge using 


H[ n+1 ] = —iV^K <3? n dt A dxi A • • • A dx n , 


( 2 . 20 ) 


where the barred coordinates are defined as follows, 


1 _ _ P~? 


t = t A|A _ 8 A * 16 , 


5 —p p—7 


X{ = Xi A _ 8 A * 16 , 

* = !,■■■ ,P 

1 ~p P+1 


Xj = Xj A _ 8 A * 16 , 

j =P+ !,-• 


( 2 . 21 ) 
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The reasonableness of the definition of the conjugate potentials will be justified later. At the 
boundary, in the equilibrium state, and ^p +4 can be expressed using A +i _ t * explicitly 
according to (2.8), 




'T>,. 11 = 


Y'VAxA -' 172 

A _) U*- A h 

aw 1 /2 /i-a_^ /2 


1 


y/2fi 


0 -a p+ 4,<t>/2 I 

A_ 


1 — Aj 


( 2 . 22 ) 


2.4 Evaluation of actions 


Since the bulk action terms and the Gibbons-Hawking term in the total action (2.1) are 
the common parts that appear in every ensemble, we define their sum as 


h = ^EH + I(f> + Ip + Ip +4 + IgH 

for later convenience. Using the solution (2.8) one can evaluate the actions, 


(2.23) 


Ic = - 


/51 / p+4^4 —p _3— 


2ft 2 


Pb 


(S-pdfi)’ +(3 -p)(A] — 2(4 — p) 


1b,p 


_3—p 

^ Pb 

/5kp+4^4—p -3—p 
2^2 Pb 

4^4—p 3—p 


(4 — p) f) -(4-P) 


A_ 


(3 -p) 


(A+A_) 

A* 


VA+ 

-5, 

-<»-*>(£ 


V+4 = - "'^2—- Pi'" < 3 - P)( A +A_)i - (3 - p) ( 


a 


= I3V P Q P $ P , 
= /5V p+ 4 Q p+ 


where 


s = lA±l5t £^ (1 _ A + )i ( 1 - = 


A 

A_ 


1 

3—p 


i-£ 

A* 


(2.25) 


is the entropy of the brane system and 

4p +4 = VpV 4 , Vp = Qx\x\ ( Pb ) ■ ■ ■ 9xpx p {Pb) Vp • Vp = j dx\ ■ ■ ■ dxp. (2.26) 

With the above results and Table 2, we obtain the classical actions for various ensembles: 
• Both charges are in canonical ensembles (we will use CC ensemble to denote this one) 


Icc — Ic T Ib,p T Ib,p +4 
/3kp+4^4 —p _3—p 

2^2 p b 
= PE-S, 


A_ 


+ (3-p)(A+A_)2 ( l + j- ) +2p-8 


-5 

(2.27) 
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where 


E = — 


kp+4^4—p _3 —p 

2^2 p b 



+ (3-p)(A + A_)2 


( i+ £) +2p - 8 


(2.28) 


is the internal energy of the brane system. This action has already been obtained in 
[30], and we include it here for completeness. 


• Dp in grand canonical ensemble & D(p + 4) in canonical ensemble (denoted as GC 
ensemble ) 


I GC =f)E-S- PV P Q P %. 


(2.29) 


• Dp in canonical ensemble &; D(p + 4) in grand canonical ensemble (denoted as CG 
ensemble) 

Icg = 0E — S — ^4p-|_4Qp-|_4^p_|_4. (2.30) 

• Both charges are in grand canonical ensemble (denoted as GG ensemble) 


Igg — @E — S — (3V p Q p <&p — pVp +i Q p+ ^p +i . (2-31) 


Here, and $ p +4 are the potentials imposed on the boundary and we stick to the con¬ 
vention that the barred quantities are the ones on the boundary. Only at equilibrium, 
Tp = 4>p and ^p +4 = <&p_|_ 4 , and 


Igc — Ic + P b,p+ 4) Icg — Ic T lb , 


.p i 


= Ic 


(2.32) 


It is easy to check that, for pb —> oo, (2.28) reduces to 


P 6 ->°o 


Vp V ^ k 2 4 P ( 4 “ P^P+ P + ( 2 “ P)P- P - (3 - p)k 


(2.33) 


the right hand side of which is exactly the ADM mass of the branes. We know that the 
free energy equals to its internal energy minus the temperature times the entropy, where 
the internal energy is the ADM mass for a gravitational system according to Bardeen et. 
al. [34]. So (2.27) and (2.33) means that we are indeed dealing with a system in canonical 
ensemble. 

Now we justify that (2.27) (2.29) (2.30) and (2.31) are indeed the correct forms of 
free energy or grand potential, and the conjugate potentials defined in (2.20) and (2.22) 
are consistent. For this purpose, we will check that the equilibrium corresponds to the 
stationary point in various ensembles. First we take the derivative of I cc with respect to 
p + and set it to zero, 


dice 

dp+ 


2vr(3 - p)Vp+ 4 n 4 _pp^_ P p b 
k 2 


/(p+i Qp; Qp+ 4) b ^(p+, Q P i Qp+ 4) 


= 0 , 


(2.34) 
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where 


b 


b(p+i Qp, Qp+ 4 ) 
/ [p+t Q P , Qp+ 4) 


Anpb' 


1 

3 — p 


(1-A+) 


1 

2 




2A+ /2 Ai /2 (l - A h 


(3 - p) A_(A_ - A + ) + 2(A+ + A^ 


+(3 -p) 


A_(A_ - A+)(A+ + A_ - 2A+A_) 

A*A+ + A*A_ - 2A + A_ 


2A + A_) 


(2.35) 


The solution to (2.34) is the stationary point of the action, which is supposed to correspond 
to the equilibrium state, and is thus expected to give the equation of state (2.19). Since 
A + < A_ < A* < 1, which can be seen from (2.9), the function / defined in the last 
equation of (2.35) is positive definite. This means the expression in the square brackets of 
(2.34) must vanish, which recovers the equation of state (2.19) exactly. This justifies our 
claim that I cc is the correct form of free energy. 

Next we take the derivative of I GC with respect to Q p and again set it to zero, 


9Igc 

dQp 


PVp 




®-®{p+,Q P ,Q P + 4 ) 


1 + 


l Qp, Qp+ 4)^ 


A*(A - A| 


I A*A + + A*A_ - 2A+A, 


PV p 


y / 2 Ke a p i)i/2 . 


$ -$(p + ,Qp,Qp+ 4) 


= 0 


(2.36) 


where 


4> = V2 ne a ^ /2 ^ p , 

In getting the first equality in (2.36), we have used the expression for Q p in (2.14), and in 
the second equality, we have used the fact that b = b(p+, Q p , Q P +i) in equilibrium state. 
The last equality in (2.36) just recovers the first equation in (2.22) where it is merely a 
definition. This proves the validity of that definition. Similarly, we calculate the derivative 


- A 


-A, 


1/2 


(2.37) 
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of I CG with respect to Q p+ 4, 


die 


PV, 


dQp+A y/ 2 ne a p + 4< ^/ 2 

2 1 


) <p - <p{p+, Qp+a) 

A — A + 


1 - 1 - 


b(p+, Q p , Q P + 4 ) \ 

* J 


+ 


3 — p A_ A*A + + A*A_ — 2A + A_ 

pv P 

y/2Ke^^/ 2 - 


p-ip(p+, Qp+i) 


(2.38) 


where 


<p = V2ne ap+4 i>p+ 4 , 


<p{p+> Q P + 4 ) = f 


Ah 

A_ 


1/2 


1 - A_ 

1 - A, 


1/2 


(2.39) 


Again in getting this result we have used (2.14) and the equation of state. Setting this 
derivative to zero would give exactly the second expression in (2.22). In the same fashion, 
if we partially differentiate I GG with respect to p + , Q p and Q p+ 4 , we would obtain all the 
three equations in (2.19) and (2.22). 

For future simplifications in the computation, we define some reduced quantities in the 
following. The reduced action in the CC ensemble is defined as 


Inn — 


k 2 Ic 


2ttpI p V p+ /fr A - p 


= -b 


A, 


1/2 


2 (^) +(3-p)(A + A_ 


> 1/2 ( 1 + i ) 


+ 2p —8 


-a - a + )‘/ 2 ( 1 - f± 


1 

3 -P 


a*; 


(2.40) 


Other reduced quantities can be defined as follows, 

, 3 ~ p 
Pb 


x=[^\ < 1 , 


Q=Q P < 1, 


_ ^ Qp +4 ^ 


3 -p 


Pb ) 


< X. 


(2.41) 


Accordingly, in the reduced variables, A +) _ ; * and the other reduced actions can be ex¬ 
pressed as 

A + = l — x, 


A_ = 1 - 


x 


A* = 


A + + A_ + y (A_ — A + ) + 4 Q A + A_ 

2(1 -Q 2 ) ' 


(2.42) 
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and 


Igc — Icc ~ (3 — p)bQ<&, 

Icg = Icc ~~ (3 — p)bq(p , 

Igg = Icc ~ (3 - p)bQ $ - (3 - p)bq<p. (2.43) 


Notice that A + only depends on x, and A_ depends on both q and x, while A* depends 
on all three variables x , q and Q. With (2.40) and (2.43), we find again the equations of 
equilibrium state corresponding to (2.34), (2.36) and (2.38), in reduced quantities, 


dice 

dQ 

dice 

dq 

dice 

dx 


-(3 -p)b{$-$(x,Q,q) 


1 + 1 - 


b(x,Q,q ) 


~(3-p)bUp-<p(x,q) 
f(x,Q,q)[b-b(x,Q,q)\, 


1 - 1 - 


b(x,Q,q ) 


\ A»(A_-A + ) 1) 

/ A*A + + A*A_ - 2A + A_J J ’ 

2 1 A_ - A + 

-^—l- -——— 

3-pA_ A*A + + A*A^-2A 


where the functions /, 6, <f> and (p have been defined in (2.35), (2.37) and (2.39). Notice 
that from (2.9), we find 0 < + < 1 and 0 < p < 1 for nonzero p- and (p + 4)-brane charges 
when the horizon is not coincident with the boundary. 


3 Thermodynamics in different ensembles 

Note that when <5 = 0, the whole system becomes a D(p+4)-brane system, and the various 
phase structures of this system have already been thoroughly analyzed in [5, 35]. Thus, in 
the following calculations, we will always assume that Q > 0 (4> > 0). Nevertheless, we 
will still compare our results with the Q = 0 case for consistency check or revealing the 
different traits in the presence of Dp-branes. 

3.1 Overview of the behaviors of b(x) 

Before performing concrete analyses in specific ensembles, we first examine a few typical 
behaviors of the function b(x, Q, q), which are useful in our later exploration. 

In the following discussions of the phase structures in these ensembles, the key problems 
to be solved are finding out the stationary points of thermodynamic potentials and which 
one is stable. That is, we need to solve the following equations, 

b{x, Q, q) = b, 

$(x,Q,q) = 

p{x,q) = ip, (3.1) 

and then determine whether the solution(s) to these equations is the minimum point of 
those thermodynamic potentials. In CC ensemble, we only need the first equation. In 



(2.44) 
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CG ensemble, we need the first and the third ones and in GC ensemble, we need the first 
and the second ones. In GG ensemble, we need all the three equations. For simplicity, 
in this paper we constraint ourselves only to the stability of branes under the change of 
the horizon size as in [6, 13, 35]. As will be shown later, this is equivalent to the thermal 
stability condition that the specific heat capacity is positive for the equilibrium state. We 
will also give a brief discussion in a later section on the other electrical stability conditions 
like those used in [7] and leave the details in a paper to appear [31]. So here we will assume 
that the second and third equations in (3.1) have already been satisfied[13, 35], and Q and 
q can be solved and be substituted in the first equation. Thus, the only problem left is to 
solve the first equation where b is a function of just one variable x, that is, 


b(x) 


= b. 


<]>=cl> or/and <p=ip 


(3.2) 


for corresponding ensembles. In that case, the local minimum condition of the thermody¬ 
namical potential would simply be 


d 2 I 
dx 2 


0 . 


(3.3) 


In equation (3.2), b is the inverse temperature on the boundary which is assumed to 
be a constant input parameter. Hence finding out the solution is equivalent to finding out 
the intersection points of curve b = b{x) and the horizontal line b = b. On the other hand, 
we will demonstrate that 

d 2 I _db 
dx 2 dx ’ 



which means the local minimum point is the intersection point where ^ (the slope) is 
negative. 

We will see in the GG ensemble, that there are only two kinds of curves for b(x) which 
are listed in Figure 1. One can see that, in the first graph of Figure 1, although there is 

b b 




Figure 1 . Typical behaviors of b{x) in GG ensembles as well as in some cases in GC and CG 
ensemble. 
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an intersection point of the two curves, it is not a minimum of grand potential since the 
slope at that point is positive. We therefore conclude that there is no stable black brane 
phase for any b given that b(x) increases monotonically. In this case, since the hot flat 
space can have the same boundary condition with the brane, the only stable phase should 
be the hot flat space. However, in the second graph, for b within some range, there could 
be two intersection points, where the slope is negative at the one with larger x. Thus we 
conclude that there will be a locally stable black brane phase for b within that range, i.e., 
the brane configuration with larger horizon is locally stable. To find out whether it is a 
globally stable phase, we have to compare the grand potential at this point with the one 
for the hot flat space, i.e. the one at x = 0. If the grand potential at this point is a global 
minimum it would be the real stable phase, otherwise it is merely a locally stable one. In 
this section, we only concentrate on locally stable phases. 

In the GC ensemble, there will be much richer phase structures. For p = 2 (D2- 
D6 system), the possible curves of b[x) are the same as the ones for the GG ensemble, 
and therefore we do not need to reanalyse them. For p = 1 case (D1-D5 system), there 
will be one more possible shape besides the two appearing in Figure 1. This new curve 
monotonically decreases as shown in Figure 2. Thus, if the constant b line intersects with 


b 



Figure 2. A typical curve in GC ensemble for p = 1 besides those in Figure 1. 

b(x) curve at some point, this point must correspond to a globally stable black brane phase. 
If p = 0 (D0-D4 system), there are three cases as shown in Figure 3. The third one as 
shown in Figure 3(c) is similar to Figure 2 with decreasing b(x). In Figure 3(a), if b = b\ 
(b > boci^max)), there is only one intersection point at which the state is stable, whereas 
if b = 62 (b < bQc(x ma x)), there will be two intersection points and only the one with 
smaller x is locally stable. For the second graph, there are more possibilities. We can see 
that if b = b\ > b max or b = 63 < b m i n , where b m i n and bmax are the local minimum and 
local maximum of b(x) as shown in the graph, there is only one intersection point and it 
is a stable black brane phase. Yet for b = 62 with b m i n < 62 < b max , there can be three 
intersection points. Denoting these three points as xi, X 2 and X 3 with x\ < X 2 < £ 3 , one 
can see that the point with x = X 2 is apparently unstable for its positive slope while the 
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Figure 3. Typical b(x) behaviors in GC ensemble for p = 0 

other two are both locally stable. As regards these two locally stable phases, the one with 
higher free energy will eventually transit to the other phase. It can be shown that there 
exists some value of b between b m i n and b max , at which the two locally stable phases have 
equal thermodynamic potentials but different entropies, which indicates a first order phase 
transition, i.e. the van der Waals-like phase transition. The second graph in Figure 3 is 
just a typical case of this kind which is similar to the one-charge black brane. In two-charge 
case b(x max ) does not reach zero, so there can be possibilities that b(x rnax ) > b m i n . We 
will discuss this case in section 3.3. 

In the CG ensemble, all three possible b(x) curves in Figure 3 also appear here. As 
argued in [30], in the CC ensemble, the delocalized Dp-branes and the D(p + 4)-branes are 
equipotent in the sense that exchanging the two kinds of brane charges would yield similar 
phase structures. This reveals some symmetry under exchange of the boundary conditions 
that we have imposed on these two kinds of branes. The GC and CG ensemble are just 
related by an exchange of the boundary conditions on these two kinds of branes and we 
will see later that the phase structures of these two ensembles are really related in this way. 

Now having found all the patterns of the b(x) curves and the locally stable phase in 
each case qualitatively, we need to be specific in each ensemble to find out in what ranges 
of the parameter Q (or <h) and q (or (p) one can have a certain pattern of b(x). Then, we 
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would know the possible state of the system given any pair of Q (or 4>) and q (or ip). The 
following subsections are devoted to this problem. 

The CC ensemble (especially when p = 1) has already been explored by Lu et al. in 
[30], and we have gathered some results in Appendix A including some facts not given 
in [30]. Interested readers are refered to their paper for more details. In the following 
subsections we will focus on the other three kinds of ensembles. 


3.2 GG ensemble 


As stated in the previous section, we adopt the one variable analysis as an embryo attempt 
to find out the thermodynamical structure. We need to solve Q and q in terms of x and the 
corresponding electric potentials on the boundary (<f> and p), from the two electromagnetic 
equilibrium equations, i.e. the second and the third equations in (3.1), using (2.37) and 
(2.39). Since the Q and q dependence comes from A_ and A*, it is convenient to first 
express A_ and A* in terms of x, d* and p, 


A_ = 


A± 

£ 


A* = A + 


1 - <h 2 


(3.5) 


where £ = ! — (! — p 2 )x < 1. Then the charges can easily be obtained, 


xp 



Q 


^(i-0 
(1 - $ 2 )^/2 ' 


(3.6) 


The reduced action which is proportional to the grand potential can then be obtained in 
terms of x, p and 


I gg = 2(4 — p)(l - y/i)[b - b 0 (x)] 


(3.7) 


where 

1 -4- Pi ( T. \ 2 1-p 

*»<*> s A (it**) < 3 - 8 ) 

We have mentioned in (2.9) that p > p 3 p > k. This inequality can be rewritten in 
terms of A+, A_ and A*, 

A + < A_ , A_ < A*. (3.9) 


The first relation in (3.9) is guaranteed by the first equality of (3.5). The second relation 
holds only for 


So, we have a restriction on x, 


x < 


1 - <h 2 

1 — <p 2 


0 < x < x 


max 


min 


i, 


1 - <h 2 1 

1 - J ' 


(3.10) 


(3.11) 
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For <& < <p, x max = 1 which means the horizon should be inside the boundary. For 
& > (p, the condition is the requirement of Q < 1. This limit should not be reached since 
otherwise A* would blow up and the size of the time direction and some space dimensions 
would shrink to zero. In fact, before it shrinks to string scale, the quantum effects should 
be large, and the supergravity approximation is invalid. Thus, if at x max the system has 
lower grand potential than at the local minimum, we will regard the system as unstable, 
since either the horizon tends to the boundary or the supergravity approximation is not 
applicable. At x = 0, the charges tends to zero, the system reduce to the “hot flat space”. 

Next, we need to find out the stationary point by setting 


0 = 


dI GG {x, <f>, p) 


dx 


where 


bccix) = 


1 




x£ 


= fGG(x) [b-b GG (x)] , 


(3.12) 


3 — p \ 1 — d> 2 


1/2 


1 - 




1-p 
2(3— p) 


Igg(x) = ( 4 -p)^7 I/| > 0. 


(3.13) 


Equation (3.12) again reduces to the familiar one, 


b = b GG (x). 


(3.14) 


Notice that at x = 0, b GG is always zero. So b GG should be increasing near x = 0 due to 
the fact that b GG > 0, i.e. db GG [x)/dx > 0 at x = 0. Finally the local stability condition 
is effectively 

d 2 I GG {x ) db GG {x) 


dx 2 dx > 

where x is the solution to equation (3.14). The derivative of b GG can be evaluated, 

A°= , bcG s r2-(5-p)(l-yVl- 

dx 2(3 — p)x £ 

Thus, the condition (3.15) is equivalent to 

2 — (5 — p)(l — <p 2 )x < 0, 


(3.15) 


(3.16) 


(3.17) 


which is exactly the same condition as for black (p + 4)-brane without p-brane charge, i.e. 
the (35) in [35], except that now the solution x depends on a new parameter 4*. 

From (3.16), we see that there is a turning point 

2 1 

■^o y ^2 

5 — p 1 — p* 

where b GG (x ) increases for x < xq and decreases for x > xq. Therefore, if x max < xq, 
b G G(,x) will be a monotonically increasing function for x € (0,x max ), otherwise it will have 
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a maximum at x = xq. For the former case (referred to as case A), the curve of boG looks 
like the one in the first graph of Figure 1, while for the latter case (referred to as case B ) 
it looks like the curve in the second graph, just as we stated in the previous section. 

We then need to figure out the condition for <1 and (p in each case. We start by 
requiring x 0 > x max , he., 

2 f 1 -l> 2 1 , 

( 5-p)(i-^) >min i 1 'T^J- (3 ' 18) 

which would give 

4> vfvf or (3 - 19> 

That means that for $ and (p satisfying (3.19), case A is applied and the system cannot 
have a stable black brane phase. On the contrary, if 


“ d ’ (320) 
bcG^x) w iH look like the curve in Figure 4. One can see that if b± > b > b 2 , there can be a 


b 



Figure 4. b G G for x 0 < x max . For $<<(>, x max = 1, otherwise x max 


i -$ 2 

i-ip 2 ■ 


locally stable black brane phase and its horizon radius is between xq and x max , where &i 
and b 2 are defined by 


bi 


b 2 


bcc{x 0 ) = 


bGGiXmax) 


2 3 -p (5 — p) 


— ^(3 pT 

(3—p)\/1—i 2 

1— p 

$(l_$2)2(3-p) 

, (3-p)y/l-if 2 


for <4> < 
for <p < <f>. 


(3.21) 


For b not in this range, there will be no stable black brane phase. 
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Combining the arguments above, we would have the diagram shown in Figure 5 which 
gives a more explicit view of the answer to the previous question. In this figure, the two 
constants are 

= ‘fo = \ [-= P ■ (3.22) 

V 5 — p 


<t> 





Figure 5. Parameter plane in GG ensemble 


Since the charges are not conserved in grand canonical ensemble, the hot flat space can 
be a possible phase. In region A, at very low temperatures with b > where dlcc/dx > 0, 
the grand potential is a monotonically increasing function of x as shown in Figure 6(a) (we 
will use the subfigure labels to denote different cases). Since the hot flat space has x = 0, 
the final stable phase should be the hot flat space. If b < b 2 , there can be two cases which 
are shown in (b) and (c) in Figure 6. In the first case (b), which has lower temperature 
than (c), the minimum of the action is still at x = 0 which corresponds to the hot flat 
space. For case (c) which has higher temperature, the global minimum of the Iqq is at the 
boundary x max which means that the hot flat space is not a global minimum of the grand 
potential and hence is unstable, and the horizon of the black brane tends to expand to the 
boundary. This means that there is no stable phase in this region. We can find out the 
condition for this case by solving the inequality Iccixmax) < 0 which gives 


^ ^unstable < 


(1—<3> 2 ) a ~P 
2(4 -p) 


l+g> 


(l-$)(l-0 2 ) 


1+ip 


for > </? 
for (p > <f> 


(3.23) 


2(4 -p) V (i^#)(iM> 2 f| 

So for temperatures larger than l/6 U nstabie the system is unstable. We do not know what 
really happens at this high temperature. Below this temperature, only hot flat space is 
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the stable phase. Region B, where max(3>,<^) < y §E^i has more cases to be considered. 
First, for very low temperatures, similar to case A, when there is no solution for (3.12), Iqq 
behaves as Figure 6 (a) and the stable phase is the hot flat space. For higher temperatures 
when (3.12) has two solutions, i.e. 62 < 6 < &i, Iqg has two stationary points, one unstable 
and the other locally stable. The possible behaviors of Iqq are shown in (d) and (e) in 
Figure 6 . For case (d) which has lower temperature, the hot flat space at x = 0 has the 
lowest grand potential, and hence is the stable phase. At higher temperatures such that the 
system corresponds to graph (e), the grand potential at the locally stable point is negative 
and becomes the global minimum. So in this case the final stable phase is the black brane. 
For certain <h and (p, case (e) can not happen. We can find out the condition for (e) to 
happen by looking at Igg{%) = 0. We put the detailed analysis in Appendix B and only 
state the results here. Firstly, when 

max($, <p) > |—— , (3.24) 

5 — p 

for 6 > 62 , Iqg(x ) is always positive which corresponds to graph (d). Under this circum¬ 
stance, the cases corresponding to graph (e) does not happen. With such $ and <p, when 
6 € (^unstable; 62 ), there is only one solution to (3.12), which also corresponds to (b) in 
Figure 6 . Hot flat space is then the global minimum as in case A. For higher temperatures 
such that 6 < fr u ns ta bi e , graph (c) also appears and the system is unstable. Secondly, when 

max(i>,<^) < ^—— , (3.25) 

5 — p 


there is a temperature To € ( 1 / 61 , 1 / 62 ) at which the locally stable point is at xo(ip) = 
pry • If the temperature is lower than To such that at the locally stable point 
x < xo(<p), Igg is positive (case (d)) and the hot flat space is the globally stable phase. 
Only for temperature higher than To, when the locally stable point x > xq (ip), the black 
brane has a negative grand potential, which corresponds to graph (e) in Figure 6 , and black 
brane is the final globally stable phase. For much higher temperature such that 6 < 62 , 
the case corresponding to graph (c) also happens, and the system is again unstable. So 
combined with case A, for all max(l, ip) > |^, the black brane can not be the final stable 
phase and for low temperatures the hot flat space is the global minimum. Only when 
max(4>,y>) < the black brane can be the final stable phase for certain temperature 


T £ (T 0 , 1 / 62 ). 

As we stated in the previous section, from the conditions (3.19, 3.20, 3.21, 3.22), we 
see that the phase structure is symmetric under the exchange of the boundary conditions 
imposed on the two kinds of branes, i.e. exchanging 4> and <p. We will also find similar 
results in GC and CG ensemble. 
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Iaa(X) 



Figure 6. Typical Igg(%) behaviors in region A and B. 

3.3 GC ensemble 

In this ensemble, the D(p + 4)-brane charge q is fixed, therefore we only need to use the 
equation 


$0,Q) = $ 


(3.26) 


to get rid of the unfixed quantity Q. Since <h is fixed at the boundary, we can solve A* 
from this equation and then Q using (2.16), 


, A + A_(l — <b 2 ) 

A + - A^ 2 ’ 

(3.27) 

d>(A_ - A + ) 

^ GC (1 - $ 2 )(A + A_) 1 /2' 

(3.28) 
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Using (3.27) and also requiring A + < A_ < A*, we can find the domain of x, 

0 < q < X < Xmax < 1, (3.29) 


where 


%max 


1 - $ 2 + <f> 2 ) 2 +4g 2 $ 2 

2 


(3.30) 


As in GG ensemble, when x — » x maxi A* tends to infinity and could not be reached in 
supergravity approximation. At x = q. at which A + = A_ = A* and Qcc = 0, the 
system reduces to extremal (p + 4)-brane. However, the horizon of the extremal brane is 
singular and the quantum gravity effect must be non-negligible, and hence the Euclidean 
action method is not applicable. So if the global minimum of the reduced action is at the 
boundary x = q or x max , we do not know what the final stable state of the system is. So, we 
regard the system as unstable in our semi-classical description when the thermodynamical 
potential at x = q or x = x max is lower than the local minimum value. 

With relation (3.28), we can find the stationary point by solving equation 


° = ~ix~ = ~ b ° c ^] ’ 


(3.31) 


where 


bcc(x) = b(x,Q GC (x,q,$),q) = 


3 — p 


1 - 


A± 

A_ 


P-1 

2(3— p) 


I A+(l-A+) 
A_(l-$ 2 ) 


r (3 -p)A_(A_ - A+) + (5-p)(A+ + A_ -2A+A_) 

JGc{x) = --—-> 0 , 


2 A+ /2 A!. /2 (1 - A+) 


(3.32) 


which gives the equilibrium condition b = bcc( x ) as before. Assuming that x = x solves 
(3.31), the stability condition becomes 


db GC {x) 


dx 


< 0 . 


Now we compute the derivative of b GG , 


db GG 


dx 


<?A 


b G c 

2x 


1 + 


(A + + A_ - 2A + A_)[2A + - (3-p)A_] 
(3-p)A+A_(A_- A+) 


(3.33) 


(3.34) 


We can see that in this expression, in the large parentheses, A* completely disappears, 
which means the equation is independent of 4*. For 4> = 0, i.e. D (j> + 4)-brane without 
Dp-brane charges, which is equivalent to setting A* = A_ in b GG {x ), (3.34) would exactly 
recover the same equation for D(p + 4)-brane in canonical ensemble ((46) in [5]). This is 
also true even for ^ 0 since it is independent of 4>, that is, (3.34) still gives the same 


condition to determine the signature of 


dbcc 

dx 


for $ = 0, which was already given as in 
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(47) of [5]. This means that the stationary point for bcc{ x ) is independent of <F and only 
depends on the charge q. To be specific, we rewrite (3.34) in terms of x and q. 


where 


db G c 

dx 


~9{x,q) 


bcc 


x 4 (l-x)(l-M(l 



g{x,q) = 


5 — p 


x 4 - 1 + 


7-p 


q x 


3 — 3p 2 2 
-r —q X 


+r 2 -p+ 


9 — 3 p 


x-(3 -p)q A . 


(3.35) 


(3.36) 


Now (3.36) is exactly the function on the left hand side of (47) in [5] (there, a different 
parameter, the “co-dimension” d which is equal to 3 — p in our case, is used). Therefore, 
condition (3.33) becomes a simpler one, 


g(x,q) > 0 . 


(3.37) 


With (3.31) and (3.37) we can proceed with the phase structure analysis for different p. 

First for p = 2 (D2-D6 system), notice that at x = q, bcc = 0 which is the same as bcG 
in GG ensemble at x = 0 . As we have mentioned in section 3.1, there are only two possible 
shapes of bcc( x ) which are similar to the cases in the GG ensemble (see Figure 1) with 
the only difference that the left end point of the bcc( x ) is at x = q here. Now we will find 
out the conditions for these two cases. From Figure 1 one can see that the difference of the 
two patterns is the sign of the slope of the curves at x rnax . So one can conclude that there 
must be a transition line on the <h — q plane on which the boc( x ) curve has dboc( x )/dx = 0 
at X max ■ Actually this is a universal feature for all p, so instead of solving it here for only 
p = 2 case we have done it for all p and put the tedious analytic computations in Appendix 
C. Here we just give the diagram in Figure 7 for p = 2, which has similar meaning as 
Figure 5. In region A, bcc is a monotonically increasing function and there is no locally 
stable black brane phase. In region B, the bcc( x ) curve is like Figure 4 and for 62 < b < b\, 
there are two intersection points and the larger one is a locally stable black brane phase. 
The boundary between A and B on q — 4> plane has the following form, 


1 - 4> 2 


Q = = 


<F(3 — 54r 


2(l-$ 2 )(3*l 2 -l), -L<$<-L. 


(3.38) 


The locally stable state may not be a global minimum of the reduced action because at 
the boundary x = q the reduced action may have smaller value. To find out the condition 
that the local minimum can not be the global minimum for 62 < b < b± one can look at 
the critical case by requiring that at b = bcc( x max ) the Icc( x max) = Igc{ x = q)- This 
equation can be solved analytically to give 


Q = 


2(1 - <b) 2 (34> - 1) 
l>(3 - 5<l) 


1 Jk 1 

3 2 


(3.39) 
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Figure 7. q — <f> parameter plane in D2-D6 GC ensemble. Region B includes all the region on the 
left of the solid curve. 

which is denoted as the dashed curve in Figure 7. On the right of this curve, when 
bo < b < b\, we always have IcG(x max ) > Igg(x) > Igg(q), hence the minimum of 
reduced action is at x = q. When b < 62 the global minimum is either at x = q or 
x = x max . There is a critical temperature To = 1/bo > 1/&2 where the reduced action at 
x = q and x = x rnax are equal. Below this temperature the global minimum is at x = q 
while above it is at x = x max ■ In both cases, we regard the system as unstable. On the 
left of the dashed curve, whether the locally stable state becomes the global minimum 
of the reduced action depends on the temperature. There is a temperature denoted as 
Testable = 1/^unstabie with 62 < Unstable < &i such that the reduced action at the locally 
stable point equals the one at x = q, i.e. Igc(x) = Igc(q)- Below this temperature, i.e. 
b > 6 uns table) the global minimum is still at x = q. Only for 62 < b < ^unstable 1 the local 
stability becomes a global one. For higher temperatures with 6 < 62, the global minimum 
of the reduced action is at x = x m ax and the system is unstable again. 

Next for p = 1, at the left end x = q, bcc( x = q) has a non-vanishing finite limit which 
is different from GG ensemble and p = 2 case. Now, g(x, q ) can be factorized as 

g{x , q) = [x + q)(x -q)(x- X-)(x - x + ), 

where 

1 + 3g 2 ± ^(l-q 2 )(l-9q 2 ) 

x ± = - - -. 

One can see that for q > 1/3, which is region C in Figure 9, and x+ are complex 
conjugates and g(x,q) > 0, which means bcc decreases monotonically for x E ( q,x max ) 
as shown in Figure 8(c). In this case, there is a globally stable black brane phase for 


(3.40) 

(3.41) 
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Figure 8. Typical bcc{x ) behaviors for D1-D5 system, (a) q < 1/3 and x max < x + , corresponds 
to region A; (b) q < 1/3 and x ma x > x+, corresponds to region B; (c )q > 1/3, corresponds to 
Region C. 


&2 < b < b\. For <7 < 1/3, it is easy to verify X- < q < x + , and therefore g(x, q ) is negative 
(bcc increases) for x € (q,x+) and is positive ( b G c decreases) for x > x + . This indicates 
that if x + > x max , bcc is monotonically increasing (Figure 8 (a)) and there is no stable 
black brane phase. This case is denoted as region A in Figure 9. For x + < x max , the 
behavior of bcc is shown in 8 (b) and there exists a locally stable black brane phase for 
62 < b < b m ax, which is denoted as region B in Figure 9. On the boundary between A and 
B, q and <f> are related by 


Q = 


yj (2 - $ 2 )( 2$ 2 - 1 ) 

3l> 


-p < d> < 1. 

V2 


(3.42) 


In region B, like in p = 2 case, the locally stable black brane state may not be the global 
minimum and the x = q point may have lower reduced action. The condition that the 
system is just to have only x = q as the global minimum of the reduced action for all 
b > 62 is 

g = I( 5 _ 2 $-|), ^<$<1, (3.43) 


which is also denoted as the dashed curve in Figure 9. On the right of the curve, the 
global minimum is either at x = q for b > 62 or at x = x max for b < 62 and the system is 
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q 



Figure 9. Parameter plane in D1-D5 GC ensemble. Region B includes all the region on the left of 
the solid line between A and B. The bcc behaviors of A, B, C correspond to (a), (b), (c) in Figure 

8. 


regarded as unstable. On the left of the dashed curve, there is a possibility that at some 
temperature T unsta bi e = 1/^unstabie the reduced action at x = q equals the one at the locally 
stable point (see Appendix E for further details on how to find this temperature). Below 
this temperature (b > ^unstable), the reduced action at x = q is the global minimum and 
the system tends to x = q similar to the cases in region C and the system is also unstable. 
The locally stable one becomes a globally stable one only when the temperature is above 
^unstable, he. b 2 < b < ^unstable- When the temperature is higher than 1 /b 2 , x = x max 
becomes the global minimum of the reduced action, and so, the system becomes unstable 
again. 

Finally, for the p = 0 case, there are three kinds of bcc{x) behaviors which are shown 
in Figure 3. Similar to the one-charge black (p + 4)-brane case [5], there exists a critical 
charge q c beyond which bcc { x ) is a monotonically decreasing function for arbitrary <f> as 
shown in Figure 3(c). When q = q c and 0 < <f> < &max, where max is obtained in 
Appendix C to be $ mai = 0.871417, bcc has an inflection point at which a second order 
phase transition occurs. At <& ma x, the position of the critical point x is equal to x rnax . For 
<f> > &max at q = q c , the critical point position x is larger than x ma x and can not be reached 
by the supergravity solution and hence bcc is also a monotonically decreasing function. 
So, in q — 4? plane, the region denoted as C in Figure 10, which contains only decreasing 
bcc(x), includes the q > q c region as well as a small region with q < q c where the solutions 
for dbac{x)/dx = 0 are larger than x m a X - In this figure, q c ~ 0.141626 is the same critical 
value as the d = 3 case in [5] as expected. So, in region C, there is always a stable black 
brane phase for b > bGc(x max ). In region A and B, curves of bcc have the same shape 
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Figure 10. Parameter plane in D0-D4 GC ensemble. Region B includes all the region on the left of 
the solid line between A and B. The bac behaviors of A, B, C correspond to (a), (b), (c) in Figure 

3. 

as in the other two diagrams (a) and (b) in Figure 3, respectively. On the boundary lines 
between A and B and between A and C, q and has following relation, 

^ ^-Ti. 2 ) ^ 3 -* 2 ^ 5 * 2 - 3 )’ Vf <4<L <344) 

In every region (A, B, C) there are chances for this system to have a locally stable black 
brane phase. In all these three cases, as x —>• q, dbcc/dx < 0 which means that the 
locally stable black brane can exist for arbitrarily small temperatures. But the minimum 
of bac does not tend to zero which means that at high enough temperature the locally 
stable black brane can not exist. In region C, at temperatures higher than 1 /bcci^max), 
the global minimum is at x max . In region A, there is a temperature T unsta ble higher than 
1 / (xmax ) i at which the reduced action at x max is equal to the one at the local minimum. 
When the temperature is higher than T unsta bie the global minimum is at x max . Under 
both these two circumstances the system is unstable and we do not know what the final 
state is. So, in region A and region C, at temperature below 1 /bac( X max ) and T unsta bi e , 
respectively, the black brane is the globally stable phase. In region B which corresponds 
to Figure 3(b), there can also be a first order phase transition under some circumstances. 
As stated in section 3.1, we need to discuss two cases separately, b m i n > bcci^max) and 
bmin < bccixmax) which correspond to the left and right regions of the dotted line in 
region B, respectively. When b m i n > bGcfamax), the system has one locally stable black 
brane phase for bQc{x max ) < b < b m i n or b > b max , which is also the global minimum. 
For b m i n < b < b max , there can be two locally stable black brane phases and the final 
phase should be the one with lower reduced action. There is a first-order phase transition 
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temperature T t = 1/bt at which the reduced actions of the two phases are equal. Above (or 
below) this temperature the larger (or the smaller) one is the global minimum. This is the 
same as the corresponding case in canonical ensemble of one-charge black brane. When 
bmin < bcc ( X max ) , we still define a temperature T t = 1/bt, at which the reduced action at 
the two local minima are equal even when one minimum has x > x max ■ In this case, we need 
to consider two subcases. The first one has b t > bGc(x max ). Then when b < bGc(x max ), 
i.e. the temperature higher than l/bGc(xmax)i the minimum of the reduced action is at 
x = x max and the system is unstable. Below this temperature, the analysis is the same 
as the above case and the system can have a first order transition at T) between the two 
locally stable minima. The second subcase is when b t < bGc{x max )- In this subcase, at 
bt the larger solution for bt = bcc(x) goes beyond x max and there is another temperature 
^unstable = 1/frunstabie(> T t ) at which iGc(xmax) equals the one at the local minimum with 
smaller x. Then, condition for x max to be the global minimum of the reduced action is 
b < ^unstable- So for temperature higher than T unsta bie the system is unstable and we can 
not say much about the phase of the system. For bGc{x max ) > b > ^unstable, the smaller 
locally stable black brane has the lower reduced action than at x = x max , and hence is 
the global minimum. So for temperature lower than T unsta bi e , the smaller locally stable 
black brane is the global minimum. The condition for bt = bGc( x max) can only be solved 
numerically and is denoted as the dashed line in Figure 10. On the dashed line, ^unstable 
coincides with bt and bGc(x max )- The left region of the dashed line corresponds to the 
bt > bGc( x max ) case which has the van der Waals-like phase structure and the right of the 
line in B corresponds to the bt < bGc(,x max ) case. 

3.4 CG ensemble 

In this ensemble, we fix the Dp-brane charge Q and the potential (p for D(p + 4) at the 
boundary. Using the “electromagnetic” equilibrium equation, 

ip(x,q) = (p, (3.45) 

the reduced action of the system depends on the only variable x, 

Icg(x) = I CG (x , Q = const., q = q CG (x , <p)), (3.46) 

where q CG is the solution to (3.45), 

q CG (x,<p) = ^. (3.47) 

This relation also guarantees that x > q CG is satisfied as long as (p < 1. At x —>• 0, qcG —> 0 
and £ —>• 1. From (3.5), we can see that A_ —>• A + , which means the curvature singularity 
is coming closer to the horizon. Actually it can be shown that the scalar curvature blows up 
at the horizon x = 0 by explicit calculations. The quantum effect near horizon, therefore, 
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may be large such that the semi-classical approach can not be applied. So, as in the GC 
ensemble, if at x = 0 or x = 1 , the system has lower reduced action than the local minimum, 
we will consider it as unstable in the semi-classical approach. 

Now we can perform the one-variable analysis by requiring 

~ = fcc(x) [b - b C c(x )] = 0, (3.48) 

where 


bcc(x ) = b(x,Q,qcc ), 

f 1 1 - A- 

JCG{x) = 


2AY 2 A^ 2 (1 - A^ 


5 - P + (3 - p)- 


A-(A - A+) 


> 0. (3.49) 


A*A + + A*A_ — 2A + A_ 

Let us assume x to be the solution of equation (3.48), then the minimum condition of free 
energy is again reduced to 

db C c(x) 


dx 

The left hand side of the above inequality is 
db C c(x ) b C G 


< 0 . 


dx 


2(3 — p)x 


2 - (3 ~A A - + 

Ar 


(3-p)A_(A^-A + ) 

A*A + + A*A_ - 2A+A_ 


(3.50) 


(3.51) 


One can easily check that when Q = 0, i.e. A* = A_, (3.51) automatically falls back to 
equation (34) in [35]. Now with (3.51) and the first equation in (3.49), we can redo all 
the analyses done in the two previous subsections. However, since the computations are 
complicated and tedious, we put the details into Appendix D for interested readers. Here 
we just list the final results. 

For p = 2 case, similar to Figure 7 in GC ensemble, we have the following graph 
(Figure 11). In area A, since bcc increases monotonically as the left graph in Figure 1, 
there is no stable black brane phase, while in area B, bcc behaves as in Figure 4, there 
exists a locally stable black brane state for 62 < b < b\. The curve separating region A and 
B is described by 


Q = 


(-1 + <^ 2 ) v / 2(! - p 2 ){2>p 2 - 1) 


x/3 <¥,< V? 


(3.52) 


ip(5tp 2 - 3) 

Again, the local minimum may or may not be the global minimum. It should compete 
with the boundary point x = 0 or x = 1. Similar to the discussion in GC ensemble, in 
Figure 11, on the right of the dashed line, the local minimum can not be the global minimum 
of the system and on the left, it can be the global minimum for certain temperature 
T G (T uns tabie, I/& 2 ), where T uns t a bie is the temperature at which the reduced action at the 
local minimum equals the one at x = 0. All the discussions in the GC ensemble can be 
used here and we would not repeat them. The dashed line is described by 


2(1 — <p 2 )(3<p 2 — 1) 11 

~ <£(3-5 p 2 ) ’ 3 < lf>< 2' 


(3.53) 
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Figure 11. Q — ip parameter plane in D2-D6 CG ensemble 

Notice if one makes the exchanges $ -H- (p and q -H- Q, (3.52) and (3.53) exchange with 
(3.38) and (3.39). So, Figure 11 is the same as Figure 7 except the labels. 

For p = 1 case, the behaviors of bcd x ) are the same as the ones in GC ensemble and 
Figure 8 can also be used here with the only difference that the domain is now 0 < x < 1. 
In the Q — (p plane, the three cases of (a), (b), (c) in Figure 8 correspond to region A, B, 
C in Figure 12. In the lower right region A, bcG is nronotonically increasing and therefore 


Q 



Figure 12. Parameter plane in D1-D5 CG ensemble 

there is no stable black brane phase. In the upper half region C, since bcd x ) is always 
decreasing, there is always a locally stable black brane state for &cg(0) > b > bed 1) 
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when constant b has an intersection point with bcG curve. In the lower left region B, for 
bcc( 1 ) < b < bmax-, where b max is the maximum of bcc , there will be a locally stable black 
brane phase since the constant b line intersects with bcc at some point where dbcc/dx < 0 . 
The curve separating region A and region B is solved in the appendix to be 


Q 


3 (p 


T2 <f<1 ' 


(3.54) 


For the local minimum in B to be a global minimum, the system should be in the 
region on the left of the dashed line. On the right of the dashed line, the global minimum 
is either at x = 0 or at x = 1 and hence the system is unstable. On the left of the dashed 
line, there is a ^unstable between bed 1 ) and the maximum b max of bcc , where the local 
minimum of the reduced action equals the reduced action at x = 0. Only for temperatures 
such that bcci. 1) < b < 6 uns tabie; the local minimum becomes the global minimum. The 
dashed line is described by 


Q = \(5-2tp-*), \<<p< 1, (3.55) 

o Lp 2 j 

One may have already found that (3.54) and (3.55) are the same as (3.42) and (3.43) except 
the exchanges ofgfiQ and <f> -H- (p. All the discussions of different phases are the same 
as the ones in GC ensemble for p = 1 and we will not repeat them here. 

The Q — <p plane graph for p = 0 case is shown in Figure 13. In all three regions, 


Q 



Figure 13. Parameter plane in D0-D4 CG ensemble 

since as x —>• 0 , bcc tends to positive infinity, there is always a stable black brane phase 
at arbitrarily low temperature. As long as (p is not zero, bcc has a finite minimum, which 
means for high enough temperatures, the black brane can not exist. To be specific, in 
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region C which occupies almost all valid area of the parameter plane, the behavior of bcc 
is decreasing similar to Figure 3(c) with the domain of x replaced by 0 < x < 1. So there 
is always a locally stable black brane phase for b > bcciX)- In region A, which takes only 
a very small portion of the plane, bcG behaves like Figure 3(a) which decreases first and 
then increases. Thus there could be a locally stable black brane phase if that solution x 
is in the decreasing segment. In region B, the bcc curve looks like Figure 3(b), and as 
discussed in the GC ensemble, there can be a van der Waals-like phase transition on the 
left of the dashed line in Figure 13 when the phase transition temperature T) = 1/bt is 
lower than 1 /&cg(1 )> where bt can be obtained by requiring the reduced actions to be equal 
at the two local minimum. On the right side of the dashed line, similar to what we did 
in the GC ensemble, we define a temperature T unsta bi e (> T t ) at which the reduced action 
at the local minimum with smaller x equals the one at x = 1. In this case, below the 
temperature T unsta bie) the smaller black brane is the global minimum. For temperature 
higher than T uns t a bie the global minimum tends to x = 1 which means the instability of the 
system and a failure of our method. A second order phase transition can happen on the 
boundary line between regions B and C, when Q = Q c = 0.141626 and (p < 0.871417. On 
this line, bcG has an inflection point where a second order phase transition occurs and the 
two phases are indistinguishable. All the detailed discussions are the same as the p = 0 
case in GC ensemble. 

As we have mentioned before, the three diagrams in this subsection are essentially 
the same as those in the previous subsection except that we substitute 4> and q for <p 
and Q. Even the critical charges are precisely equal, q c = Q c . This indeed supports the 
statement that the Dp branes and the D(p + 4) branes are equipotent as far as only the 
thermodynamics is concerned, which has been pointed out in [30]. One may think of it as 
a symmetry, 


tp -H- <h, Q «->• q. 


(3.56) 


4 More general thermodynamic stability conditions 


The stability condition (3.3) for the thermodynamic potential we used here is the one with 
respect to the fluctuation of the horizon size. It can be shown that this condition is the same 
as the positivity of the specific heat capacity for each ensemble as follows: According to the 
first thermodynamic law and our definitions of the potentials, for equilibrium state cLE = 
TdS + (3-p)$dQ + (3-p)pdq , one can find (ff ) Qq = T(x,Q,q ) (§§) Qq with T(x,Q,q ) = 
1 /b(x,Q,q), which can also be checked by explicit calculation. For GC ensemble, from 
Ice = bE(x,Q,q ) — S{x,Q,q) and by using the equilibrium condition 


0 = 


dice 

dx 


Qq 







(4.1) 
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Comparing with (2.44), one finds that ± (§§) Qq = l (tffOg,, K f( x ’Q’V) > °- The 
specific heat for CC ensemble is defined as CQ q = T(| t)qq = — So f° r the 
specific capacity heat to be positive, (J|)q < 0 should be satisfied which is the result of 
the stability condition. Similarly, for the other ensembles, we have 


f dI G c \ 

\ dx ) 


1 fdS\ 


Qv 


= (*-%*: =-(b-b)c Q9 


b \dx J 


Q<p 


db\ 

9x )q <p " 


= 0 




( s j2i\ = ( 5 - 6)4 


\ & /$ 




<l>g 

S).r*- 


.0 


Qq 


(CG); 

(GC); 

(GG), (4.2) 


where the equilibrium conditions for <p and in (3.1) are used. Comparing with (3.31), 
(3.48), (3.12), and using (3.49), (3.32), and (3.13), one can find the equivalence of the 
positive specific heat capacity conditions and our previous used conditions (3.33), (3.50) 
and (3.15) for different ensembles, respectively. 

The positivity of the specific heat describes the stability of the system under the 
fluctuations of horizon size or, equivalently, the tenrperature/entropy. In [4, 7] the electric 
fluctuations are also considered and the stability condition of one-charge black holes to 
these fluctuations depends on the positivity of the electric permittivity. Similarly, here 
since the Dp-D(p + 4) system has two kinds of charges, there can be two electrical stability 
conditions for this system when these two kinds of electric fluctuations are considered. 
According to the thermodynamics, the general stability conditions could be deduced either 
from the maximal entropy or minimal energy criterion. As a result, one can find out the 
requirements for different “response functions” suitable for the discussion of stabilities in 
various ensembles, 



£t<p = 


^Sq = 


£t<p = 


£sq = 



> 0 

> 0 

> 0 

> 0 


(CC); 

(GC); 

(CG); 

(GG). (4.3) 


If there is only Dp brane charge, the positivity of £ is not needed. In [35], in the grand 
canonical ensemble for one-charge black brane, the electric stability condition is consistent 
with the thermal stability condition of positive specific heat. However, this is no longer true 
in our case if one of the Dp- or D(p + 4)-branes (or both of them) is in canonical ensemble. 
The electrical stability analysis will reshape the phase structures in GC/CG/CC ensembles 
and the branch of the smaller black brane will be unstable in all possible cases when the 
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van der Waals-like phase transition is supposed to happen. This is similar to the result in 
[7], where part of the smaller AdS black hole branch including the van der Waals phase 
transition point is unstable when electrical stability condition is considered. After all, the 
brane system is not so simple as the liquid-gas system because carrying charges is the 
nature of branes and thus the electrical stability is almost inevitably as important as its 
thermal stability. The complete analysis on the electric stability of this system is a little 
complicated and we shall defer the details to a companion paper [31]. 

5 Conclusions and outlook 

In this paper, we have discussed different thermodynamical ensembles of Dp-D(p + 4) 
system, where p = 0,1,2. The two kinds of charges can be in either canonical or grand 
canonical ensemble separately, so there can be CC, GG, GC, CG ensembles. CC ensemble 
has already been discussed in [30] and we focus on the other three ensembles in this paper. 
In the following we summarize the results obtained in this paper. 

For GG ensemble, the potentials for D p and D(p + 4) are fixed. At very low tempera¬ 
tures, the hot flat space is the stable phase. Depending on the values of and <p, at higher 
temperatures black brane phase may or may not be a globally stable phase. For larger $ 
or <p which satisfies (3.24), black brane can not be a globally stable phase. In this case, 
below T unsta bie the globally stable phase is the hot flat space and above it the horizon of the 
system approaches the boundary and the quantum effect will be important, therefore we do 
not know what happens in this system under these circumstances. In this case, the system 
is regarded as unstable in our semi-classical approach. Only for small <f> and ip satisfying 
(3.25), the black brane can be a globally stable phase at temperature T E (To, l/b(xm ax )). 
At much higher temperatures, like in case A, the horizon of the system also tends to the 
boundary and the black brane is unstable. As in most grand canonical systems, there is 
no van der Waals-like phase transition. 

For GC ensemble, in which the D(p + 4) charge q and D p potential are fixed, the 
D2-D6, D1-D5 and D0-D4 behave differently. In D2-D6 system, on the right side of the 
dashed line in Figure 7, the global minimum of the reduced action is either at x = q or at 
x = x max , and the black brane can not be a stable phase. On the left of the dashed line, in 
a comparably small region of 4>, the black brane phase can be the final stable phase only 
when the temperature is in range (T unsta bi e , 1 /bcci^max))- Below or above this range, the 
global minimum is at x = q or x = x max , respectively, and the system is thus unstable. 
For D1-D5 system, in the region below q c = 1/3 in Figure 7, the discussion is similar to 
the D2-D6 system. For the region above q c , the black brane is always the final stable phase 
for temperature lower than 1 /bGc{x max ) and higher than 1/bcciQ)- For higher or lower 
temperatures, the global minimum for the reduced action is at x = x max or x = q which 
means the instability of the system. For D0-D4 system, there can be a van der Waals-like 
first order phase transition between a small black brane and a larger one in the region on 
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the left of the dashed line in region B in Figure 10, which is below the critical charge q c 
at which a second order phase transition happens. The critical charge is independent of <h 
and is the same as the one for the black D4-brane in the canonical ensemble. Below the 
first order phase transition temperature T t (or above T t and below 1 /b(x max )), the final 
system is in a small black brane phase (or a large one). For higher temperatures, the global 
minimum is at x = x max again, which means that the system is unstable and is beyond 
our approach to handle. In the region on the right of the dashed line in B, or in region A, 
for temperatures below T uns tabie or 1 /b(x max ), respectively, the black brane is always the 
final stable phase, though there may be a larger metastable one in the corresponding part 
in region B for some temperature T € {l/b max , 1 /b(x max )). In region C, for temperatures 
lower than 1 /b(x max ), the final stable phase is always the black brane. 

The fact that the critical charge q c in D0-D4 GC ensemble is independent of <F is an 
unexpected result. Notice that according to (3.28), Qcc still depends on 4* non-trivially. 
Recall that for extremal Dp-D(p + 4) branes satisfying the Harmonic function rules and 
preserving 1/4 supersymmetries, there is no binding energy [36], which means that the two 
kinds of charges do not affect each other. However, we are considering non-extrenral cases 
here, and the two kinds of charges must be correlated, which is already demonstrated in 
the CC ensemble where the critial line correlates both charges as in Figure 14. Whether 
this independence of of the critical charge in GC ensemble means that something special 
happens in such a critial condition is an interesting problem for future work. 

From above discussion, we can see that D2-D6 GC ensemble behaves more like GG 
ensemble. The bcc and boo have the same behavior as described by Figure 1 and changes 
similarly as we tune the <f> in both GC ensemble and GG ensemble at small (p . The 
difference is that at x = 0 in grand canonical ensemble there is the hot flat space but at 
x = q in GC ensemble the extremal black brane has naked singularity and our method 
fails. There could be new phases emerging around x = q in GC ensemble like hot flat space 
in GG ensemble due to quantum effects. D0-D4 GC ensemble is more like the canonical 
ensemble. There can be a van der Waals-like phase transition and there is a critical charge. 
Tuning q is similar to that in the one-charge black brane canonical ensemble for small "F. 
But for large the van der Waals-like phase transition may disappear. 

In this paper, we only consider the thermal stability condition, i.e. the positivity of the 
specific heat. The phase structure for CG ensemble is almost the same as the one for GC 
ensemble except for the interchange in q «->• Q and $ <p. Similar to CC ensemble [30], 

the phase structure of GG ensemble is already symmetric itself under this interchange. So 
the smeared Dp charges and the D(p + 4) charges in the Dp-D(p + 4) system are equipotent 
at least so far as only the thermal stability is concerned. We will show explicitly in another 
paper [31] that the electrical stability condition will still preserve this symmetry, though 
more constraints shall be imposed. However, we will show there that the phase structure 
will also be modified just like the situation which has been handled in [7] for charged AdS 
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Figure 14. Critical lines in CC ensemble. The left diagram corresponds to D0-D4 case and the 
right one corresponds to D1-D5 case. The right diagram has already been shown in [30]. Both 
critical lines are computed numerically. 

black holes. 


Appendices 

A Phase structures in CC ensemble 

Compared with the other ensembles, the analysis in CC ensemble is much simpler due to 
the fact that both charges Q and q are fixed, which leads to the consequence that 

bcc( 1) = 0. (A.l) 

Hence the shape of Figure 3 (a) would never occur, i.e., systems that have van der Waals-like 
structures can only have bcc( x ) with shapes shown in Figure 3 (b) and (c). Therefore, there 
is only one transition line in the Q — q plane which is also the critical line on which a second 
order phase transition occurs. However, unlike in GC or CG ensemble where only D0-D4 
can have van der Waals-like phase transitions, in CC ensemble, both D0-D4 and D1-D5 can 
have this behavior. In Figure 14 we show the critical lines in Q — q parameter plane. In both 
diagrams, B corresponds to the region where bcc( x ) is a monotonically decreasing function 
and A corresponds to the region where van der Waals-like phase transition can happen. 
In D0-D4 case, the critical line intersects with the coordinate axes at Q c = q c = 0.141626, 
which exactly matches the result in the presence of only Dp or D(p + 4) charges [5]. 

For D2-D6 system, there are no critical behaviors. However, depending on where ( Q , q) 
pair lies in the Q — q plane there still can be two possible shapes of bcc( x )> as shown in 
Figure 15. The diagram for Q—q plane is shown in Figure 16, in which region A corresponds 
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Figure 15. Possible shapes of bcc(%) for D2-D6-brane in CC ensemble. 

q 



Figure 16. Q — q plane for D2-D6-brane in CC ensemble. 


to the left diagram in Figure 15 and region B corresponds to the right one. The boundary 
line between these two regions can be obtained analytically to be 


Q 

q 4Q-1 


or 


Q 


Q 

4q-V 


(A.2) 


where ^ < Q,q < 1, which shows the symmetry between charge Q and q. 


B The condition for globally stable black brane phases in the GG en¬ 
semble 


When x = x we have 

J 5 -p( x \ 1/2 (1-0^ A-i/2 3 -p\ 

GG 3-p\i-&) 1 +e 1 / 2 v 5 -p)' 


(B.l) 
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where £ = 1 — x(l — ip 2 ). Now we need to find out when I GG < 0, or when 

2 


£< 


3 — p 
5 — p 


This leads to the condition that 


4(4 — p\ l 

x > x 0 {(p) = ——or b < b GG {x 0 (<p)). 


(B.2) 


(B.3) 


(5 — p) 2 1 — ip 2 

However, we have a restriction for x which is xo < x < x max . That means only if xq{ 0) < 
Xmax we can have x > x${(p). It is easy to check that xq(( p) > xo is automatically satisfied, 
whereas for xo((p) < x rnax to hold we need 

3 ~P 


max{<h,99} < 

5 — p 

Thus we find the condition for the black brane phase to be globally stable, 

3 — p 


max{^>, (p} < 


5 — p 


and b < b GG {x Q ((p)). 


(B.4) 


(B.5) 


C Parameter planes in GC ensemble 

We have argued in section 3.3 that we can find the transition lines between region A and 
region B in Figure 7 for p = 2 by setting the b GG curve at the right end point to be flat. 
This can also be used in arbitrary p. That is, 


diXmaxi Q) — 0 , 

where x max satisfies following equation (see (3.30)), 

X 2 max ~ I 1 - ® 2 )Xmax ~ q ^ = 0. 


(C.l) 


(C.2) 


The above quadratic equation has two solutions and x max is the larger one. Therefore, 
there is another restriction on x m(7T , 


Xmax ^ 


1 - & 


From (C.2), we can express q through x r 


^2 _ ^max 


- (1 - $ 2 )x r , 


4)2 


(C.3) 


(C.4) 


Putting this relation into (C.l), we can factorize g(x max ,q) as follows, 

X max(l — Xmax ) 


9^x m ax i q) — 


2 <h 4 


-{ [3(3 — p) - (7 -p)<h 2 ] x r 


-2(3 -p - <h 2 )(l - 


(C.5) 
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Thus we find the relation between x max and <f> when (C.l) is satisfied, 


_ 2(3-p-$ 2 )(l-$ 2 ) 

Xmax ~ 3(3-p)-(7-p)4 > 2 ' 

Combining (C.4) and (C. 6 ), we finally reaches the relation between q and 4>, 

1 - 4> 2 


(C. 6 ) 


Q=sT\ 


2(3 — p — 4> 2 )[(5 — p)4 > 2 — (3 — p)]. 


$|3(3-p) - (7-p)$ 2 | 

From (C. 6 ) and (C.3), the restriction on x max becomes a restriction on 4>, 


(C.7) 


4 > < 


' 3(3 -p) 
7 — p 


(C. 8 ) 


For p = 0 or 1, this restriction is satisfied automatically since $ < 1 while for p = 2, this 
inequality reduces to $ < yj 3/5. Now it is easy to see that, for p = 0 or 1, (C.7) reduces 
to (3.44) or (3.42) and for p = 2 together with (C. 8 ) this recovers (3.38). 

For p = 1, there is another transition line between region B and region C in Figure 9 
which appears when the hoc curve becomes flat at the left end point. That is equivalent 
to setting 


dbgc(q) 

dx 


(g ~ V3) 2 

q 2 { i - q) 2 


boc(q) = o. 


(C.9) 


This gives a condition which is independent of 4), 


<7 = 1/3. (C.10) 

So the transition line is the horizontal line shown in Figure 9. 

For p = 0, there is a critical line between region A and C on which the bcc h as a 
inflection point and the first and second order derivatives of bcc vanishes at the same 
point x. That means, we only need to combine two equations, 


g(x,q) = 0 , 
dg{x,g) = 
dx 

This set of equations are solved by 

x = 0.292675, q c = 0.141626. 


(C.ll) 


(C.12) 


Since x < x max , this critical point exists only when x max > x, which means there is an 
upper bound for <3>, 


Tmax = 0.871417. 


(C.13) 










D Parameter planes in CG ensemble 


We have argued that the behavior of bcc{x ) depends on the signature of b' CG {x). Hence we 
need to find the zero points of b GG (x). With (3.47), one can easily find following relations, 



(D.l) 


where 


_ i+e+v(i -zy+*qh 

2(1 -Q 2 ) 


With this relation 


one can evaluate b GG (x) and set it to zero, 



bcc Z [3(3 - p) ~ (7 - p)£] - A(£ + 1)(3 - p - 2£) 
2(3 — p)x £(£A + A-2£) 


which results in, 


. £[3(3-p)-(7-p)£] 

(£ + l)(3-p-2£) ' 

The second definition in (D.2) is equivalent to requiring 

0 = (1 — Q 2 ) A 2 — (1 + £)A + £> and A > ^ 


by the same token of (C.2) and (C.3). Combining (D.4) and (D.5), one gets 


(D.2) 


(D.3) 


(D.4) 


(D.5) 


_ l~g / 2(3-p-£)[(5-p)£-(3-p)]" 

' 3(3-p)-(7-p)eV e 

< f 1, for p = 0,1 
' \ f, for p = 2 


(D.6) 


By the same spirit of the previous appendix section, we examine when the slope at 
each end point of b G c curve changes its signature. That is to set x = 0 and x = 1 for left 
and right end points respectively. For x = 0 to be the stationary point, or equivalently at 
£ = 1, which only makes sense when p = 1 (because only this case we have finite limit of 
bcc)-, we have Q = 1/3. This corresponds to the horizontal line in Figure 12. For x = 1, 
or equivalently £ = <p 2 , we get 


Q = 


l -^ 2 


3(3 - p) - (7 - p)<p 




2(3 -p- y 2 )[(5 -p)ip 2 - (3 -p)] 
V 2 


(D.7) 
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and (p < y | for p = 2. This relation represents all curves in Figure 11, 12 and 13. 

Finally, we calculate the critical charge Q c for the p = 0 case. From (D.6) we know 

that 


Q c = 


!-<£ /2(3-0(5£-3) 


9 - 7<e V £ 

Another critical condition is that b GG (x) = 0. Since 


b cc( x ) 


(4£ - 1)A - (14£ - 9) + 


(e + i)(2e-3)(A-i) 




we find another equation for Q c , 

2(1 _ g 2 ) = « + im- 3 ) 


1 + 


£(7£ - 9) 

Combining (D.8) and (D.10), we get only one sensible solution 

£ c =* 0.759367, Q c = 0.141626. 


2(l-Q2)A-l-e 

(g-l)(2g-3)(5£-3) ~ 

- 3)(n<e - 9) 


(D.8) 


(D.9) 


(D.10) 


(DU) 


This also says that, on the critical line, x and <p are related through a simple formula, 

1 - £ c 0.240633 _ , 

x = --§ = --—• (D.12) 

1 — <p z 1 — </r 

Since x < 1, this formula also indicates that (p has a maximal value beyond which Q c does 
not exist, i.e., (p m ax = \[Q. — 0.871417. 


E b 


'unstable 


in GC ensemble and the dashed lines in Figure 7 and 9 


Our starting point is to find the critical b = 6 uns tabie above which we have Igc(q ) < Igc( x ) 
and the system at x = x is not globally stable. As a bonus of this process we will find the 
dashed lines in Figure 7 and 9. When b = & uns tabie we have 


Igc(q) — Igc( x ) 

where x takes the value such that 

b = b{x). 

In (E.l) the left and right hand side can be evaluated, 
Icc{q) = (3 -p)bq, 

Igc(x) = b 8 - 2p — (5 — p) ^ 

A + (x) 


(E.l) 

(E.2) 


8-2p- (5-p) f xww l - ( 3 ~P) (A+(a;)A_(x)) 1/2 


1 - 4* 2 


1/2/ A ^\3/2 


= b 


8 — 2p — 2 


1 - 


A+(x) 


A_( 


A_(x) 


- (3 -p) [A+(x)A_(®)] 2 - (3 -p) 


A(x) 

1 i 

2 

.A+(®). 



(E.3) 
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where we have used (2.43) (3.27) (3.28) and (E.2) to obtain the final expressions. Combining 
(E.l) and (E.3) we get an equation of x, 


8 - 2p - 2 


A+(x) 

A_(x) 


- (3 -p) [A+ (x) A_ (x)] 2 - (3 - p) 


A-(x) 

_ A +(x) 


= (3 -p)q. (E.4) 


That is, given any pair of (&,q) we can solve (at least numerically) the above equation to 
find the x that makes I G c(q) and Igc(x) equal. One can easily see that x = q is always a 
solution to this equation which is not the solution we want. However, it can be shown that 
if there exists a solution other than q, this solution is unique and therefore there exists a 
unique 6 uns tabie = b(x). Consequently, there may be regions and boundaries thereof in the 
g-d* plane within which this solution always exists. These regions should lie in region B of 
Figure 7 and 9 and the boundaries are exactly the dashed lines in these figures. We have 
argued that there is an upper limit for x which is x max defined in (3.30), therefore on the 
boundary we have x = x max . Bringing this specific value of x in (E.4), we find 


= 2 (l-$)(3-p-$)[(5-p)$-(3-p)] 

Q 3 -p T[3(3-p)-(7-p)T] ' 1 ' ; 

Setting p = 2 and p = 1 respectively the expected functions in (3.39) and (3.43) could be 
recovered. 
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